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Abstract 

We present a web interface which allows to conveniently setup calculations based on the 
BioFET-SIM model. With the interface, the signal of a BioFET sensor can be calculated 
depending on its parameters, as well as the signal dependence on pH. As an illustration, 
two case studies are presented. In the first case, a generic peptide with opposite charges 
on both ends is inverted in orientation on a semiconducting nanowire surface leading to a 
corresponding change in sign of the computed sensitivity of the device. In the second case, 
the binding of an antibody/ antigen complex on the nanowire surface is studied in terms of 
orientation and analyte/nanowire surface distance. We demonstrate how the BioFET-SIM 
web interface can aid in the understanding of experimental data and postulate alternative 
ways of antibody/antigen orientation on the nanowire surface. 



Author Summary 
Introduction 

A bionanosensor is most generally described as a device that allows the detection of an 
analyte (e.g. 11+ ions, small molecules, proteins, DNA, viruses, cells) at ambient condi- 
tions where the dimensionality of the sensitive component is on the nanometer scale. The 
sensitive component can be either a functionalized nanotube, nanoribbon or nanowire, 
the latter being the focus of this paper. Currently, a large research effort is dedicated to 
the development and application of bionanosensors including pH measurement [l] , protein 
sensing [2]-[5], DNA detection [6|[7], blood analysis [s], nanotechnology based medicine jo]. 
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and the description of fundamental performance limits of these sensors 10 12. A number 



of reviews describe the bionanosensor 13 -17 and its components. In addition to the ex- 
perimental work, simulators of bionanosensors are being developed and several numerical 



models have been presented 18 22 



Most simulators are aimed at providing a measure of the current or conduction through 
the sensor, which are the prime experimental targets. This requires, in principle, the 
description of the charge distribution on the sensor and within. From the charge dis- 
tribution, the potential within the sensor is calculated which in turn is required for the 
calculation of the current. The calculation of the potential can be either numerical or 
analytical. 

In this paper, we present a computational tool to simulate a bionanosensor which is based 



on an analytical model 23 25 and which can calculate the sensitivity of the nanosensor 
and the pH dependence of the signal upon binding of a protein. The use of an analytical 
model is mainly motivated by the fact that this model does not require extensive compu- 
tations but still allows to gain a qualitative understanding of the biosensor problem in a 
straightforward manner. 



Furthermore, we have demonstrated 24,25 that 1) the experimental data can be repro- 
duced with sufficient accuracy to help interpret them and 2) going beyond the simplifi- 
cations inherent in the model may not be warranted until the key properties of current 
BioFET experimental set-ups are known with greater precision. We note that the pre- 
sented method, which we refer to as BioFET-SIM, has gained popularity in the biosensing 



community and is being actively incorporated into present day research 26 -29 



Because of the reduced required computational effort, it is possible to incorporate the 
model into a browser based application which by doing so can be made accessible to a 
wide range of users. Our goal is to provide a tool from which indications for trends in 
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predictions can be obtained with minimum effort of preparation and time. To further 
improve the usabihty, the model is coupled to an atomic representation of the protein 
structure in a way many researchers in the biocomputational field are familiar with. Such 
an application is an ideal tool for gaining insight and obtaining semi-quantitative solu- 
tions to the problems at hand which can be of valuable guidance in the design process of 
an experiment, for optimization of experimental parameters and rationalization. 



We relate our application to other simulators where we point out the BioSensorLab 30 



which implements settling time, sensitivity and selectivity of the biosensor, Nanowire 



31 , which allows to carry out self-consistent three dimensional simulations of a silicon 



nanowire or Medici [32 j, a commercial simulator. Custom prepared simulators 18 have 
also been described. To the best of our knowledge, out of all available simulators targeted 
at modeling of biosensors, the tool we present in this paper is the first to combine a three 
dimensional visual representation of the biomolecule to be studied directly in the browser 
with a method to solve the biosensor problem. 

The primary use of our tool is to model the binding of proteins to nanowire surfaces for 
which the structure is available in the PDB database. However, using a custom prepared 
structure, it is also possible to model binding of an antibody/antigen complexes [2^, an 
illustration of which is provided in the results section. The authors further envision the 
application of the program to the modeling of DNA binding. 

To put the use of this application into context, we note that every program in general 
requires a certain amount of preparation. When using atomic detailed structures, the 
molecular structure of the pH dependent charge distribution has to be generated, which 
usually involves the combined usage of a number of different software tools, each dedicated 
to a particular task. If different orientations of the analyte charge distribution are to be 
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studied, the procedure needs to be repeated for each orientation. Furthermore, for each 
orientation, the nanowire is covered differently, the evaluation of which requires additional 
manual effort. 

The motivation for the development of the web interface is to eliminate this effort as 
far as possible. The effort of assigning partial charges to the amino acids is essentially 
removed. Instead, using the web interface, any number of orientations of analyte towards 
the nanowire surface can be generated within minutes. Furthermore, the coverage of 
the nanowire is calculated instantaneously for any orientation of the analyte towards the 
nanowire surface. Allowing the user to interactively adjust the orientation of the analyte 



through a Jmol 33 applet provides a maximum of visual feedback and allows to generate 



the coordinates of the charge distribution as straightforwardly as possible. 

In addition to the web interface, a command line version of the program is available. The 
command line version is used together with a special BioFET-SIM input file written by 
the interface which can be used to reproduce a given calculation locally. 

The BioFET-SIM Online web interface is hosted at www.biofetsim.org, the source code 
for both the web interface and the command line version of the application is hosted at 
Github, the URL to the repository is found on the interface page. 



Methods 

A BioFET nanosensor consists mainly of a semiconducting nanomaterial covered by an 
oxide layer and a (bio-)functionalization layer. The device is usually immersed in an 
electrolyte containing the analyte. In the following, we describe the implementation of 
each of these domains in the BioFET-SIM program. 
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Theoretical Background 

The sensitivity AG/Go, where AG is the difference between the conductance upon binding 
G and the base conductance Gq, of the nanowire is evaluated using a Thomas- Fermi 



screening model for the charge carriers in the nanowire 34-36 . We point out that in 



this model, only one type of carrier in the entire nanowire is considered and the nanowire 
material is assumed to resemble a low density metal. 

In this context, the description of the electrostatic problem of the biosensor is governed 
by two major assumptions regarding the carrier concentration which are 1) the carrier 
concentration is assumed to follow an uniform distribution throughout the wire and 2) it 
is assumed not to be influenced by the electrostatic potential due to surface charges. The 
sensitivity (assuming a p-type doped nanowire) is evaluated by 



AG 



Repo 



-r 



where R is the radius of the nanowire, e is the elementary charge, Pq is the hole carrier 
density in the nanowire and o";,. is the corresponding surface charge density of the charge 
i on the biomolecule b containing m ionized sites (residues and termini) [24]. li^tot is 
the distance of the discrete charge Qi above the nanowire surface, which will further be 
discussed below. Ti is given by 
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In Eq. [2| I is again the distance between the discrete charge and the nanowire surface, 
i.e li^tot in Eq. [ijand is the Debye screening length of the electrolyte/buffer solution 



(the expression for which is given below). The expression for F is found in section SI. 2 of 
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the supporting material but for the purpose of the discussion can be considered a factor 
with values ranging from zero to unity. F and F/ are dimensionless functions quantifying 
the actual sensitivity of the nanowire (F) and the effect of a^^ (F^) and arise from the 
solution to the Poisson equation in cylindrical coordinates given the boundary conditions 



of the problem 25 . F depends on both Xd (describing the ionic strength of the buffer) 
and the Thomas- Fermi screening length Xtf (describing the electric field screening within 
the wire), whereas F; depends only on Ad. 

The screening model for the wire is a simplification in the sense that possible deactivation 



of dopants at the surface 37 or the increased dopand concentration near the surface 



compared to the semiconductor bulk 38 is solely described by the screening length Xtf- 
For a p-type (n-type) semiconductor, the screening length Xtf is related to the charge 
carrier density po {tiq) through 



-^TF - J / 1/3 (3) 

where is the relative permittivity of the nanowire material and m* is the effective mass 
of the charge carrier (po would be replaced by hq for an n-type nanowire). From Eq. 
[3} we note that Xtf can be interpreted as a measure for the charge carrier density in 
the nanowire under no applied bias and that therefore this parameter can be used to 
simulate the effect of the back gate in an experimental setup. We note that the accuracy 
of the predicted signal is strongly dependent on the quality of the estimation of the charge 
carrier density in the wire, thus for best predictivity, this parameter has to be as close to 
the actual value of the experimental setup as possible |4]. 

We further note that the described linearized model is not capable of describing non-linear 
effects such as inversion mode of operation. However, the model distinguishes between 
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accumulation/depletion mode of operation by allowing to choose between a n- or p-type 
material and different values of \tf- 



The oxide layer is known from earlier studies 24 to have an important effect on the 
predicted sensitivity and is a key component of a BioFET sensor. The gate dielectric is 
understood to be in part responsible for biosensor degradation due to the incorporation 
of charges when exposed to solvent (through ion diffusion) |39|. However, in our approach 
surface charges formed on the oxide layer surface and within are not taken into account. 
In other words, only the signal generated by a charged system bound at the surface of 
the sensor is considered. The change in signal given by, e.g., a change in pH which can 
affect the surface charge density of the oxide layer, is considered as background signal. 
The biofunctionalization layer is currently implemented solely as a distance parameter, 
providing a measure of the spatial extension of the linker molecule. Charges on the linker 
molecules are not considered. By using the same distance between the surface and the 
sensed protein for all proteins, we imply that all proteins are binding in one orientation 
to the nanowire surface. This is being further discussed below. 

We note that in principle the surface functionality of the nanowire is non-uniform |2j and 
requires a combined description of the pH dependent charge on the linker molecules as 
well as the oxide where a common description of the charge of the oxide layer is through 



the site-binding model 40 



The influence of buffer characteristics on device performance has been described 41 , 42 



and we note that the electric screening of the analyte by the buffer can have a consid- 



erable effect on the predicted signal 24 . As stated above, the screening of the analyte 



signal by the electrolyte is implemented through the expression Ti which depends on the 



Debye length A/j = \/^2nX^' where ks, T and A^^ indicate, the Boltzmann constant. 
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temperature and Avogadros constant, respectively. The expression for the ionic strength 
is given by / = 1/2 Cizf where q indicates the concentration of ion species i and Zi is 
its formal charge. Furthermore, eo and denote the free space dielectric constant and 
the relative permittivity of the electrolyte, respectively. 

The description of the electrolyte by the given approach assumes 1) that the electrolyte 
is in equilibrium, i.e. the chemical potential is at a minimum and 2) that the value for Ad 
used in the expression for Ti, Eq. [2} is equal to the Debye length of the electrolyte. We 



note that in principle these values can differ due to the biofunctionalization layer 23 



The enzyme protonation states are described classically. Depending on the pK^ value, 
the charge on residue i is calculated as a function of pH using Eq. |4] 

where p{i) = 1 for i G {Asp, Glu, C-, Tyr, Cys} and p{i) = else (the charge is evaluated 
only for ionizable residues). In Eq. lil qi{pH) can be interpreted as the probability of the 



amino acid being protonated 43 . The three-dimensional protein charge distribution is 
obtained from placing the charge calculated from Eq. |4] at the average of the coordinates 
of the terminal atoms of the side chain of residue i. The charges of the enzyme residues are 
calculated solely depending on the pH of the electrolyte and their respective pK^ values 
as computed by PROPKA. Binding to the nanowire is assumed not to affect these pK^ 
values nor to disrupt the overall protein conformation. 

Interface Operation 

The interface is shown in Fig. [1} The interface operation is grouped into three steps: 
1) Initialization, 2) Jmol based calculation setup and 3) BioFET-SIM-signal/pH-response 
calculation. 
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Initialization, Fig. [T]A.. On loading the interface, the user is requested to grant access 
to the client computer by the Java applet. This is required if the user wants to be able 
to save a Jmol state file or to restore a previous session. 

The calculation is prepared by setting the PDB identifier and the pH value. Alternatively, 
the user can upload a custom made molecular structure (in PDB format), which is then 
being submitted to the web interface. In case the user uploads a custom prepared PDB 
file to the web interface, this PDB file has to contain the MODEL and END tags, a generic 



example is provided in section SI. 5 of the supplementary material. After successfully 
uploading a PDB file, the structure can be loaded into the interface by using its file name 
(without extension) in the PDB identifier input field. 

The following steps are carried out in the background by clicking "Initialize" . The server 
first checks the availability of the requested PDB file in an internal database (assuming no 
file was uploaded) and downloads the PDB file of the biological assembly from the PDB 



database |44] (www.pdb.org) if needed. The file is processed using PDB2PQR vl.7 45 ,46 



to fix any missing side chain atoms. The structure is realigned to its main rotational axes 



and its center of mass is placed at the coordinate origin using the VMD 47 packages 
ORIENT and lal . 0. The pK^ values of the ionizable amino acids are computed using 



PROPKA v3.0 48 . Since ligand molecules are discarded from the PDB file during the 
preparation of the calculation, the additional computational effort of calculating the pK^ 
values by PROPKA v3.1 |49| can be avoided. The C-terminus is added by the PDB2PQR 
routine (in form of an OXT atom), while the backbone nitrogen of the first amino acid of 
each chain represents the N-terminus. In order to display the generated discrete charge 
distribution, the charges and the respective coordinates are written to a PQR file where 
atomic radii are arbitrarily set to l.OA. This PQR file thus contains only as many entries 
as there are ionizable residues and backbone termini present in the biomolecule. After 
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carrying out these steps, the structure is loaded into the Jmol applet. The CPU time 
required to carry out all of the above described steps depends mostly on the size of the 
molecule. On average, a time of 1-2 minutes is observed for a PDB file representing 
a medium sized protein (around 300 residues). The most time demanding step is the 
realignment of the structure to the coordinate axis. However, if a PDB identifier is 
selected for which the aligned structure is already present on the server, the realignment 
step is skipped and the time requirement is significantly reduced. Using the "Reinitialize" 
checkbox, the interface can be instructed to carry out all previous steps even if a structure 
with the same name is already present on the server. This is required if a file is uploaded 
for which an older version with the same name is already present on the server. 

Calculation setup, Fig. [Tp, C. A Jmol representation of the computed charge distri- 
bution overlayed with a ribbon representation of the biomolecule is displayed. A flat plane 
of carbon atoms illustrates the nanowire surface (without having any influence on the com- 
puted results). In this Jmol applet, the user can adjust the orientation of the biomolecule 
towards the nanowire surface allowing to take into account how the biomolecule binds 
according to the position of its binding sites. Also, it is possible to study the effect of 
different orientations on the signal, in particular if a specific orientation has a significantly 
different signal compared to other orientations. The parameters (Fig. [T]C) defining the 
BioFET-SIM calculation can be adjusted below the Jmol applet and they correspond to 
the parameters introduced in Table [Tj Recommended lower and upper limits for the pa- 
rameters, as well as a tool to calculate the charge carrier density from Xtf (Eq. |3]), is 
provided on a separate help page, the link for which is found on the interface. 
The BioFET-SIM calculation requires the computation of the normal distance (z-coordinate) 
of the discrete charges from the nanowire. Since the structure is placed at the coordinate 
origin, the atoms and charges have formally positive and negative 2;-values in the Jmol 
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applet (Fig. |2]A.). When submitting the calculation by clicking the "BioFET-SIM" button, 
internally all charges are offset by the most negative z-value, Zmin, Fig. |2]B. Due to the 
offsetting of the coordinates, any free space between the biomolecule and the nanowire 
introduced by adjusting the orientation has no effect on the computed results. Together 
with the biolinker- and oxide layer thickness, the total distance, k^tot, of each discrete 
charge from the nanowire surface is computed (Fig. [2p) and used in the evaluation of 
the sensitivity by Eq. [Tj 

Upon submission of the calculation to the server, the number of biomolecules covering the 
nanowire in the given orientation is determined by dividing the nanowire surface with the 
area of the face of the bounding box of the biomolecule oriented towards the nanowire. 
Fig. |3} In doing so, it is assumed that the nanowire is completely covered by biomolecules, 
that all biomolecules are oriented in the same way and, as stated above, that all bound 
biomolecules are equally distant from the nanowire surface. This is reasonable to assume, 
when considering high affinity binding between as e.g. in the biotin and (strept-)avidin 



complexes 50 . Complete coverage of the nanowire has been demonstrated experimen- 
tally fsTj. Alternatively, the web interface also allows the manual setting of a parameter 
defining the number of proteins covering the nanowire surface independent of the ori- 
entation of the biomolecule or the nanowire surface area. This feature is added to the 
web interface because it is questionable if the number of molecules should adjust with 
orientation or not. For non-globular proteins, the required area on the surface can vary 
strongly with orientation, however the number of linker molecules is assumed to remain 
the same for two different orientations. 

By selecting the "Single" option, the web interface also allows to use the single charge 



model 24 , where the overall charge of the analyte is placed at the geometrical center of 
the enclosing bounding box and the discrete charge distribution within the protein is not 
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considered explicitely. The single charge mode of interface operation is useful when no 
particular binding orientation is favoured. 

When a calculation has been carried out, a Jmol state file can be saved on the user ma- 
chine. This file allows to restore a session at a later point in time. As stated above, this 
option is only available if the user grants access to the signed applet, else the state file 
can not be written to the user machine. We demonstrate the restoration of a session in 
an instruction video (URL is found on the interface page). 

Calculation of results, Fig. [ip. Two types of calculations can be performed: 

1. BioFET-SIM signal, giving the sensitivity as a function of the parameter selected 
using the "Plot" radio button in the indicated range (this calculation type is illus- 
trated in the discussion of the generic peptide model) 

2. pH response, giving the sensitivity as a function of pH for the parameters entered 
(shown in Fig. [T]D) 

The pH response signal is computed by evaluating the BioFET-SIM signal at different pH 
values which will correspond to different partial charges on the residues of the protein. 
The plotted data and a specially formatted input file for the command line version of 
the BioFET-SIM program can be downloaded after the calculation is carried out. The 
input file contains all parameters together with the charge distribution and allows to 
carry out the calculation with the command line version of the BioFET-SIM program 
(system requirements and usage instructions are given in supplementary material). For 
convenience, a label indicates the sensitivity and the base conductance computed at the 
given set of parameters. Fig. [Tp, top. 
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BioFET-SIM Command Line Version Description 

The command line version of the BioFET-SIM program can be used to run calculations 
locally after the orientation of a biomolecule towards the nanowire has been established us- 
ing the web interface. In order to do so, a BioFET-SIM input file (with .bfs file extension) 
containing the charge distribution and the BioFET-SIM parameters can be downloaded 
from the web interface after running a calculation. The input file is in binary format and 
not directly human readable. However, using the command line version of the BioFET- 
SIM program, the parameters can be viewed and adjusted. The command line version 
of the program can be used for automated calculations. The usage of the command line 



version is illustrated in section SI. 5 of the supporting material. The command line version 



is open source and is hosted at Github (the URL is provided on the web interface page). 



Results and Discussion 

To illustrate the use of the web interface, we perform two case studies. In the first case, 
a generic linear peptide is placed on the nanowire and the dependence of the sign of 
sensitivity on the orientation of this peptide is evaluated (Figs. |4]^-C). In the second 
case, we demonstrate the effect of different orientations of an antibody/antigen complex 



on a relative sensitivitiy value and relate to experimental work by the Reed group 27 . 



Generic Peptide Model 



The generic peptide used in this study is prepared using the molecular building feature 
of the PyMOL 52 program. The peptide consists of two (protonated) Lys at the N- 
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terminus and two (deprotonated) Asp residues at the C-terminus which are bridged by 8 
Ala residues (the termini contribute the third charge at each end of the molecule). The 
overall charge is -0.23 formal charges at pH 7.4, the nanowire configuration corresponds to 
the default values as shown in Table [ij the calculation is carried out for a p-type nanowire. 

In the orientation of Fig. |4]A. the negatively charged aspartic acids are close to the nanowire 
surface, in Fig. |4j3 the positive and negative charges are roughly equally distant from the 
surface, and in Fig. |4p the positively charged lysine residues are close to the nanowire 
surface, respectively. For each orientation, the dependence of sensitivity on Debye length 
Xd is computed and shown in Fig. [5]A.. 

It is clearly visible how the orientation affects the sign of the sensitivity. When the negative 
charges on the Asp residues are closer to the wire (Fig. |4]\), positive charge carriers are 
accumulating in the wire leading to increased conductivity. When both Asp and Lys 
residues are equally distant from the wire (Fig. |4|3), the effect on the charge carriers 
cancels. When the Lys residues are closest to the wire (Fig. |4p) the situation is reversed 
such that positive charge carriers in the nanowire are repelled by the positive charges on 
the peptide, rendering the nanowire in depletion. The slightly different absolute values of 
the sensitivity at a given value of for the two vertical orientations are due to the not 
exactly inverted orientation of the peptide on the nanowire and due to the fact that the 
charges are not distributed in a perfectly symmetrical way on the peptide (the Lys side 
chains being longer than the Asp side chains). This results in slightly different population 
numbers on the nanowire for the two orientations. 

The observed signal is further rationalized in terms of the functional form of F;, Eq. |2j 
In Fig. |5p, F; is plotted as a function of the charge-surface distance li^tot for different 
values of Xd- The plots illustrate that F^ is comprised in the [1, 0] interval. When 
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considering the orientation of the generic peptide (~ 3.8 nm long) reported in Fig. |4]\, the 
aspartate charges are close to the surface, which means Ti^Asp is close to 1 and contributes 
significantly to T through the product F;. ■ o";,. (Eq. |2]). The lysine side chain charges, 
instead, are at a distance from the surface for which Ti^iys is observed to be close to zero. 
Therefore F^. ■ at. of the lysines is minimal. Only by diluting the buffer solvent (e.g. 
Xd = 4.0 nm) these charges could contribute more to the signal. 



Antibody Study 

In the second case study, the web interface is used to study the effect of binding differ- 
ent orientations of an antibody/antigen model complex. Experimentally, it was shown 
that different orientations of the antibody are responsible for different signals, which 
are indicative of different distances between the charged antigen and the surface of the 
nanowire [27|. Two possible binding states of the antibody appear plausible. In one state, 
the antibody is bound by an N-terminus which is located on the antigen-binding fragment 
(Fab), Fig. |6} In the other state, the antibody is bound to the nanowire surface by one or 
both C-termini at the base of the antibody. When binding through the N-terminus, the 
antigen is reported to bind at a distance of 5.9±0.6 nm and when binding through the 
C-termini, the antigen is reported to bind at a distance of 8.4±0.4 nm above the nanowire 



surface 27 



Using the web interface, different orientations of the antibody/antigen complex have been 
generated for both states and studied in terms of their effect on sensitivity. A description 
of the preparation of the molecular model of the antibody/antigen complex used for the 



study and molecular images of the different orientations are provided in section SI. 3 and 



Fig. SI and the raw data is reported in Tables SI and S2 of the supplementary material 
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We note that for the purpose of this study, the antibody is considered not to interact 
electrostatically with the nanowire. It is introduced merely as an advanced form of spacer 
and as a guidance in the construction of different orientation schemes for the binding. 

The two states are characterized by different degrees of freedom to orient the anti- 
body/antigen complex on the surface. For both states, the considered orientations are 
defined in the schemes of Fig. [7} 

The orientations A-C correspond to the orientations considered in the experimental work 



by Reed et al 27 . The orientations D-G were generated while considering further possible 



binding orientations under the conditions imposed by the binding through N- or C-termini. 

In the orientations A, C and G, the C-termini at the base of the antibody restrict the 
movement of the antibody with respect to the nanowire. In contrast, when binding 
through the N-termini on the Fab, the antibody is more free in its movement on the 
nanowire surface and the antigen can be placed at a larger range of different distances. 
These orientations are indicated by the schemes B, D, E and F. 



Following the derivation provided in section Sl.l of the supporting material, the average 



distance between the antigen and the nanowire surface, /, can be estimated by fitting the 
expression for the relative sensitivity factor, F^/r™'^^, for different values of A^. For the 
orientations A-G of the antibody bound in one of the two states, the computed values of 
the relative sensitivity factor are shown in Figs. [8]A. and B together with the value of I 
obtained from the curve fit. 

Considering the limited movement of the antibody when binding by the C-termini, the 
antibody is required to remain upright on the nanowire surface. For the orientations A, C 
and G, I is found to be between 9.8 and 14.5 nm. For the state bound by the N-terminus 
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on the Fab (giving rise to the orientations B, D, E and F), / is found to be in the range 
of 5.9 to 17.8 nm. In this range, the lowest two values of /, 5.9 and 7.9 nm, correspond to 
orientations in which the antibody is lying on the surface of the nanowire (orientations 
E and F). The molecular image of the antibody orientation corresponding to curve 1^ is 
shown in Fig. |9j 

The orientation F results in a computed relative sensitivity factor for which the fitted 
value of / is in best agreement with the antigen/nanowire surface distance of around 5.9 
nm reported in the experiment. 

From this case study, it is postulated that in addition to the orientations considered so 
far by Reed et al., a number of other orientations appear plausible as well. Based on our 
findings, we postulate that orientations E or F of the antibody are most likely to explain 
a signal corresponding to an antigen/nanowire distance of 5-7 nm and an upright position 
(A, G) is most likely to explain a signal corresponding to an antigen/nanowire distance 
of 10-13 nm. In addition, we observe that the orientation C is likely to correspond to an 
orientation where the antigen is placed even further away from the nanowire surface and 
is thus unlikely to explain the experimentally observed low value of /. 

Conclusions 

We describe a web interface to model the signal of protein binding to a nanowire based 
BioFET sensor. 

In the model, the nanowire is described using Thomas-Fermi theory, assuming uniform 
carrier distribution of one carrier type and no deactivation of dopants. The oxide layer 
is described through its thickness and permittivity, without considering surface or buried 
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charges. The biofunctionaUzation layer is considered to provide a distance measure of 
the analyte to the nanowire surface, however it is not considered as carrying charges 
and is assumed to bind all analytes identically. The electrolyte is described using Debye 
theory assuming equilibrium conditions. The charge distribution on the analyte (protein) 
is calculated from PROPKA and is assumed not to be influenced by the binding to the 
nanowire surface. 

We point out that the presented method is considered a tool which can provide qualitative 
insight into the biosensor problem, especially in cases where not all key experimental 
parameters are available |24|. 

The web interface presented in this work enables efficient and convenient use of the 
BioFET-SIM model. The automated generation of the pH dependent charge distribu- 
tion and the freely rotatable 3D representation of the biomolecule allow to study the 
effect of geometrical orientation and charge distribution on the sensitivity. By providing 
these features, the web interface significantly reduces the previously required manual ef- 
fort of preparing a BioFET-SIM calculation. In addition, the web interface is platform 
independent making it possible to use the BioFET-SIM model within any operating sys- 
tem environment and requiring only a Java enabled web browser being installed on the 
local machine. 

A specially formatted input file prepared by the web interface allows to redo a calculation 
using the command line version of the BioFET-SIM program locally. 

For studying less complex systems consisting of only one formal charge (e.g. binding of 
glutamate), it is also possible to use the previous version of the web interface. 

Two applications of the web interface are illustrated. In the first, the change in sign of the 
sensitivity is demonstrated using a generic linear peptide model with opposing charges on 
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each end. In the second apphcation, the web interface is used to study the binding of an 
antibody/antigen complex. A number of orientations are studied and we use the web in- 
terface to interpret experimental data published by Reed et al [27] . Based on the findings, 
it is concluded that the previously postulated orientation of the antibody/ antigen complex 
is not necessarily the most reasonable explanation of the observed signal. It is postulated 
that an orientation where the antibody/antigen complex is lying on the nanowire surface, 
is most appropriate to explain the observed value of the antigen/nanowire distance re- 
ported by the Reed group. Furthermore, based on our findings, we rule out one of the 
proposed orientations as not plausible. 
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Figure Legends 
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Figure 1. BioFET-SIM Web Interface. A: Upload or request of protein structure 
and pH setting. B: Jmol visualization of protein on nanowire surface. C: BioFET-SIM 
parameter section. D: BioFET-SIM calculation result. 
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Figure 2. Definition of distance reference system. A: Protein center of mass 
aligned to coordinate origin (z-axis is offset to left for clarity). B: Protein structure 
offset by z-min. C: Definition of distance of discrete charge, k^tot, to NW surface. 




Figure 3. Illustration of occupied surface area on NW. Red area indicating 
coverage of the NW by a single biomolecule. 
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Figure 4. Manual rotation of a generic KK{8A}DD peptide in the Jmol 
applet. A: Asp close to NW, A=136976, B: Asp and Lys equally distant from NW, 
A=25462, C: Lys close to NW, A=139821. 
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Figure 5. Dependence of sensitivity on orientation and Debye length and F; 
dependence on k^tot- A: The data series corresponds to either the Asp or Lys residues 
being close to the nanowire surface in Figs. |4j4. and C. The black data series corresponds 
to both Asp and Lys residues being equally distant from the surface, as in Fig. |4]B. B: 
Dependence of F^ on k^tot, Eq. |2| for different values of and R = 10.0 nm. For 
Xd = 0.2 nm, the function vanishes for li^tot > 3.5 nm. 




Figure 6. Illustration of custom prepared antibody /antigen system. 

C-termini in gray at antibody base, N-terminus on Fab in brown. 
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Figure 7. Definition of studied antibody orientations. The antigen is indicated 
by a red arc. The antibody base is indicated by a double line. The point of attachement 
to the NW surface is indicated by the small circle. Orientations A, C, G are bound by 
the C-termini. Orientations B, D, E, F are bound by the N-terminus. 
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Figure 8. Computed dependence of relative sensitivity factor on Debye 
screening length. A: Binding by C-termini. B: Binding by N-terminus. 
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Figure 9. Antibody in orientation F lying on NW surface. Antigen charge 
distribution is bound to the right Fab. 
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Tables 



Domain 


Parameter 


Default 


Unit 


Description 


NW Properties 


Lnw 


2000 


nm 


JNW length 




Rnw 


10 


nm 


NW radius 




Xtf 


2.04 


nm 


Thomas-Fermi screening length 




ei 


12.0 




NW permittivity 




M 


lE-2 


m^V^is-i 


Charge carrier mobility, 










/i = 2^(37r^n)^/'^, n: electron concentration 




Kq 


1.11E24 


m~^ 


Charge carrier density, k, G {n, p\, k, ;= k{\tf) 




K 


P 




NW doping type, K e {n, p} 


Oxide layer and 




2.0 


nm 


Oxide layer thickness 


biolinker properties 


£2 


3.9 




Oxide layer permittivity 




lb 


1.0 


nm 


Biolinker thickness 


Solvent properties 


\d 


2.0 


nm 


Solvent Debye length 




£3 


78 




Solvent permittivity 


Biomolecule properties 


TV 


4000 




Number of biomolecules on NW 










(computed internally or defined by user) 



Table 1. BioFET-SIM parameters. The analytical expression for X^f is given in 
Eq.g 
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Electronic Supplementary Information 

BioFET-SIM Web Interface: Implementation and Two Applications 
Martin R. Hediger, Jan H. Jensen, Luca De Vico 

Sl.l Derivation of Relative Nanowire Sensitivity Factor 

According to our model, the sensitivity of a generic p-type nanowire, in the single charge 
approximation, can be expressed as 

AC 2 
and the base conductance Gq can be expressed as 

Go = ^^^^ (S2) 

In these equations e is the elementary charge, R is the radius of the nanowire and po the 
hole density, fi is the charge carrier mobility. F and F/ are dimensionless functions quan- 
tifying the actual sensitivity of the nanowire and they depend, among other parameters, 
on the distance / of the sensed charge from the nanowire surface and the buffer Debye 
length Ad [l] . When using a multiple charge model, we interpret / as the average distance 
of the sensed charges from the nanowire surface. 

If the physical and geometrical properties of the nanowire are fixed. Go is constant. If we 
only consider the sensing of charges immersed in the buffer (i.e. let as = 0), Eq. SI can 
be simplified and it is possible to express the change in conductivity as 

AG = KTTiab (S3) 
where K collects all constant values. Using the expression for F/ 
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we can now define a value for tlie buffer Debye length at maximum dilution, A^"^, and 
express the maximum change in conductivity at this value as 

When considering a highly diluted buffer, i.e. A^"^ ^ /, we can express Tf^"-^ as 



(S6) 



The ratio between the change in conductivity at a given Debye length and the maximum 
possible value becomes 




AG KTTiab TTi 



and after reordering it is possible to obtain 

AG r""''^ 

After inserting the explicit expressions, we obtain 



max 



(S7) 



(S8) 



2i&(l + ^exp(//A^; 



(S9) 



where we define r^/rj""^ as the relative sensitivity factor. 

Using Eq. 2 from previously published work [ij, it is possible to compute the values of F 
for different Debye lengths, including T""'^ for A^""" = lOOOnm. Using BioFET-SIM it is 
possible to obtain the value of AG (and AG™"^) simply by multiplying ^ with Gq. It 



is then possible to plot the l.h.s. of Eq. S8 for different values of \d- This plot can be 



fitted to the r.h.s. of Eq. S9 where / is the fitting paramenter. From a series of measures 
at different Debye lengths, it is then possible to obtain the average distance of the sensed 
charge from the nanowire surface. 
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SI. 2 Expression for F 

For simplification of notation, in tlie expression for F, we use the tliickness t := Rnw + ^: 
Tlien, F is given by 



t 



ei 



Xtf 



\^ Xtf 
(SIO) 



In Eq. SIO , Jq, h, Kq and Ki are the modified Bessel functions of first and second kind, 
respectively [1]. 



SI. 3 Antibody and Antigen Preparation 

A suitable complex structure of a generic antibody and an antigen used by Vacic et 
al. [2] was prepared. Only few full antibody structures have been resolved. We used 
the structure of an intact IgG2a monoclonal antibody, ascension code IIGT [3|. For the 
antigen we used the structure of the SEA domain of human mucin 1, with ascension code 
MUCl |4j. 

The antigen structure was rigidly docked to the antigen-binding site of the antibody using 
AutoDock fsj and visually checked with the program PyMOL [6] for a reasonable docking. 
The scope of this docking was only to obtain a feasible complex structure. 

Since we were interested in the sensing of only the antigen, the antibody structure had 
to be made as neutral as possible. In order to make the antibody as neutral as possible, 
all positions in the antibody sequence were mutated to glycine using PyMOL. When 
BioFET-SIM computes the charges of a protein, a positive charge is assigned to the N- 
terminus and a negative charge to the C-terminus, depending on the corresponding pK^ 
values calculated by PROPKA [?]. In order to counter balance the charges of the termini, 
in each of the four chains of the antibody we mutated the residue at the N-termini to 
aspartate and the residue at the C-termini to arginine. 



SI. 4 Antigen Sensing 

The BioFET-SIM parameters were set as in Tab. [T] With these parameters we obtain 
Go = 279.0 nS. 
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The different orientations of the neutral antibody/antigen model complex which were 
tested are shown in Fig. SI A biofunctionalization layer of 0.5 and l.Onm was added 
for C- and N-terminus binding to the nanowire surface, respectively. We considered a 
pH of 7.4 and kept the number of proteins fixed to 4000 protein units. The BioFET- 
SIM results for all orientations A-G at the values of Debye length employed by Vacic et 
(3.07, 9.7 and 30.7nm) and X'B"' 



al 



lOOOnm are reported in Tab. SI 



S2, the derived values for the l.h.s 



together with 
of Eq 



S8 



the corresponding values for F. In Tab 
are reported. The data from Tab. S2 is plotted in Fig. S2, together with the fitting 
parameter. 
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Figure SI. Different orientations of the neutral antibody/antigen complex. In 
orientations A, C and G, the complex is bound by the C-termini, in orientations B, D, E 
and F, the complex is bound by the N-terminus. 




ESI-6 



Table SI. Sensitivity 



Orientations 


A 


B 


C 


D 


E 


F 


G 




Debye length [nm] 


AG 
Gn 


AG 
Gn 


AG 
Gn 


AG 
Gn 


AG 

Gn 


AG 
Gn 


AG 

Gn 


r 


3.07 


0.005466 


0.006495 


0.019521 


0.040992 


0.120386 


0.254474 


0.056443 


0.0514 


9.7 


0.549386 


0.260311 


0.403622 


0.666037 


1.11689 


1.574997 


0.877657 


0.1257 


30.7 


1.866226 


1.161079 


1.508332 


2.083992 


2.981374 


3.83572 


2.529128 


0.2403 


1000.0 


5.24021 


3.506654 


4.356305 


5.737461 


7.834166 


9.786054 


6.795953 


0.5651 



Table S2. Relative sensitivity factor, values for orientations A-G. 



Debye length [nm] 


A 


B 


C 


D 


E 


F 


G 


3.07 
9.7 
30.7 


0.01146788 
0.471323318 
0.837504368 


0.020363313 
0.333725685 
0.778646598 


0.049265849 
0.416530283 
0.814235818 


0.078549161 
0.521877579 
0.85417693 


0.168945023 
0.640925665 
0.894942693 


0.285889363 
0.723539345 
0.921745075 


0.09131075 
0.580583226 
0.875169108 
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SI. 5 Command Line Version of BioFET-SIM Usage 

PDB Format Requirements: 

The uploaded PDB file is required to contain the MODEL, TER and END tags. In addition, 
individual chains of the structure are required to contain the chain label (i.e. A, B, ...). 

MODEL 1 

ATOM IN ASN A 2 0.209 -1.748 -0.613 1.00 0.00 

TER 

ATOM IN HIS B 2 12.057 2.821 9.469 1.00 0.00 

END 

A number of instruction videos are available, the links to these arc found on the interface 
website, www.biofetsim.org. The link to the command line version repository is found at 
the same URL. 
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System Requirements: 

The command line version of the BioFET-SIM program requires the following applications 
and libraries to be installed on the local host. 

• Python 2.5 or higher (not Python 3) 

• Numpy and Scipy libraries installed 

Basic usage is demonstrated by issuing 

[user] $ python bio_run.py 
BioFET-SIM usage: 

$ python bio_run.py — calc <input.bfs> 
or 

$ python bio_run.py — set <param> <val> <input.bfs> 



Starting a BioFET-SIM calculation using the input file "kkSadd.bfs": 

[user] $ python bio_run.py — calc kkSadd.bfs 

# BioFET-SIM Calculation 

# Date of calculation: 2012-03-21 17:17:51 

# Calculation target: kkSadd 

# pH: 7.4 

# Comment: BFS Input generated by interface. 



Adjustable Parameters: 

L_d 2.0 

L_tf 2 . 04 

eps_l 12.0 

eps_2 3.9 

eps_3 78.0 

lay.bf 1.0 

lay_ox 2 . 

mu 0.01 
n_0 l.lle+24 

num.prot 4000 

num_qi 6 

nw.len 2000.0 

nw_rad 10.0 

nw_type P 

target kkSadd 
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Base Conductance [nS] : 279.352916413 
Sensitivity: 0.115675065297 

In the output, the labels have the following meaning (with the corresponding symbol 
given m Tab. [T| Qi given in Eq. |4]): L_d: Debye screening length Xd; L_tf : Thomas- Fermi 
screening length Xtf', eps_l: Nanowire permittivity ei, eps_2: Oxide layer permittivity 
€2] eps_3: Solvent permittivity 63; lay_bf : Biolayer thickness 4; lay_ox: Oxide layer 
thickness lox] niu: Charge carrier mobility fi; n_0: Charge carrier density kq] nuin_prot: 
Number of biomolecules A^; num_qi: Number of charged sites within each biomolecule 
Qi] nw_len: Nanowire length Lj^fw] nw_rad: Nanowire radius Rnw'-, nw_type: Nanowire 
doping type K; target: PDB identifier of the studied biomolecule. 
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A parameter can be adjusted by the "-set" option, followed by the label for the corre- 
sponding parameter and the new value: 

[user] $ python bio_run.py — set L_d 3.0 kkSadd.bfs 
Parameter adjusted. 
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